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We present a search for the pair production of scalar top quarks (ii), the lightest supersymmetric 
partners of the top quarks, in pp collisions at a center-of-mass energy of 1.96 TeV, using data 
corresponding to an integrated luminosity of 7.3 fb _1 collected with the DO experiment at the 
Fermilab Tevatron Collider. Each scalar top quark is assumed to decay into a b quark, a charged 
lepton, and a scalar neutrino (p). We investigate final states arising from t\t%— > bbp,rvv and titi — > 
bbTTVV. With no significant excess of events observed above the background expected from the 
standard model, we set exclusion limits on this production process in the (m^mp) plane. 

PACS numbers: 14.80.Ly, 12.60.Jv, 13.85.Rm 



Supersymmetry (SUSY) pQ is a space-time symme- 
try that associates a bosonic partner with each standard 
model (SM) fermion, and a fermionic counterpart to each 
SM boson. The mass eigenstates of the scalar fermions, 
fi and /2 , are the results of the mixing of the SUSY part- 
ners of the chiral states Jr and fi, . The mass splitting be- 
tween fx and /a depends on the mass of the corresponding 
fermion. It has been suggested [5] that the large mass of 
the top quark (t) can induce a large splitting between the 
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two stop-mass eigenstates with the consequence that the 
lightest scalar top quark t\ may be sufficiently light to be 
produced abundantly at the Fermilab Tevatron Collider. 
If i?-parity [3] is conserved, scalar top quarks would be 
produced in pairs in pp collisions, either through gluon 
fusion or quark-antiquark annihilation. 

In the minimal supersymmetric extension of the stan- 
dard model (MSSM) [4 with i?-parity conserved, squarks 
(q) usually decay directly into q — > q Xi where the 
lightest neutralino Xi 1S the lightest supersymmetric par- 
ticle (LSP). At the Tevatron, this mode is kinemati- 
cally disfavored for the lightest scalar top quark, be- 
cause of the large mass of the top quark. If in addi- 
tion, ra; < mj+m-t , where xt is the lightest chargino, 

the decay channel t\ — » bxl is not accessible, and the 
only two-body decay that would be allowed is the fla- 
vor changing decay ii—t cxi [5]. Searches related to this 
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channel have been reported by the ALEPH, DELPHI, L3 
and OPAL Collaborations [6]. The CDF and DO 
Collaborations have searched for scalar top quarks in fi- 
nal states with acoplanar charm-jets and large imbalance 
in transverse momentum {I$t)- Other possible decays 
are the three-body modes t\ — > bWxi, *i — ^ bH + Xi, 
t\ — > b£v, and t\ — > b£v, where H + is the charged Higgs 
boson, and v and I are the sneutrinos and sleptons, su- 
perpartners of the neutrinos and leptons, respectively. A 
possible four-body decay, t\ — > , where / repre- 

sents a fermion, mediated by virtual top quark, chargino, 
sbottom, slepton and first /second generation squark ex- 
change, has also been suggested [9]. Searches related to 
this channel have been reported by the DO Collabora- 
tion [TU]. It is thought that the three-body decay modes 
may be important and even dominate the loop-induced 
cxi mode [TT]. Searches for scalar top quark pair pro- 
duction in bb££ vv final states have been reported by the 
ALEPH, L3, and OPAL Collaborations [6]. DO [10l [12j- 
IT5] and CDF [TB] have searched for scalar top quark pairs 
in the final states bb££ vv, with leptons in e/x, mi, or ee 
channels. No Tevatron searches have yet considered sig- 
natures with hadronically decaying r leptons although 
SUSY could well appear at the Tevatron in final states 
with taus [17] . 

In this letter, we search for scalar top quark pair 
production in events with r leptons, assuming that the 
branching fraction B(t\ — > b£v) = 1 and that the sneu- 
trino is either the LSP or decays invisibly. We search for 
stop pair production through the decay tiii — > bb^irvv, 
or t\ti — ¥ 66rr(— ¥ fivv)vv, in a data sample correspond- 
ing to an integrated luminosity of 7.3 fb _1 at a center-of- 
mass energy of 1.96 TeV, collected with the DO detector 
at the Fermilab Tevatron pp Collider between April 2002 
and July 2010. The signal topology consists of one iso- 
lated muon, one isolated r lepton, and I$t coming mainly 
from undetected sneutrinos, and unreconstructed or mis- 
measured jets. 
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FIG. 1: Diagram contributing to the three-body decay 
blv. 

The three-body decay t\ — » b£v proceeds mainly 
through a virtual chargino xt (Fig- [l])- In the MSSM, 
Xi is the lightest mass eigenstate of the charged gaugino- 
higgsino mass matrix that is a mixing of the wino 
W + and the higgsino H + , SUSY partners of the W bo- 
son and the charged Higgs boson, respectively. If xt 1S 
wino-like, the leptonic decay xt ~~ ^ occurs with 



equal rate to all lepton flavors. If xt is higgsino-like, 
the decay xt 1S enhanced, owing to the large 

Yukawa coupling of the r lepton. In that case, the decay 
ti — > brv can be dominant. We consider two scenar- 
ios in our search that depend on the composition of the 
chargino. The wino scenario is defined by B(t\ — > b/iv) = 
B(t\ — > brv) = 1/3. For the higgsino scenario, we choose 
B(ti — > b[iv) = 0.1 and B(i\ — > brv) — 0.8, which cor- 
respond to the maximal values reached with a scan of 
the MSSM parameter space using SUSY-hit [15]. In both 
scenarios, the signal is a combination of the bbfirvv and 
bbTTvv final states. 

The DO detector [T9T[2T] is designed to optimize detec- 
tion and identification of particles arising from pp inter- 
actions and comprises dedicated subsystems surround- 
ing the interaction point. The central tracker resides 
within a liquid-argon/uranium sampling calorimeter and 
muon detectors. Charged particles are reconstructed us- 
ing multi-layer silicon detectors and eight double layers of 
scintillating fibers in a 1.9 T magnetic field produced by 
a superconducting solenoid. After passing through the 
calorimeter, muons are identified using 1.9 T toroids and 
a muon system composed of three layers of drift tubes 
and scintillation counters. Events are selected for offline 
analysis through a three-level trigger system. All events 
contributing to this analysis are required to pass one of a 
suite of single-muon triggers based on information from 
the tracking and muon systems. 

For each event, the best primary vertex (pv) is selected 
from all the possible reconstructed interaction vertices as 
the one with smallest probability of originating from a 
minimum-bias interaction [22] . To ensure efficient recon- 
struction, the location of the primary vertex along the 
beam direction is restricted to \z pv \ < 60 cm, where z pv 
is the longitudinal position with respect to the center of 
the detector. 

Using central track segments pointing to hit patterns 
in the muon system, muons are identified in the region 
|?7| < 1.8, where -q is the pseudorapidity [23] • Their tra- 
jectories are required to have both drift-tube and scin- 
tillator hits that match a track in the central tracker. 
Muons that are not isolated are rejected if the sum of 
the transverse momenta of tracks inside a cone of ra- 
dius 1Z = yj (A0) 2 + (Ary) 2 = 0.5 around each muon (4> 
being the azimuth), divided by the transverse momen- 
tum of the muon, is less than 0.15. The sum of the 
transverse energies in the calorimeter in an annulus of 
0.1 < 1Z < 0.4 around the muon, divided by pt^,, is 
also required to be less than 0.15. Only muons with p^ > 
15 GeV are considered in the analysis. A veto on cosmic 
ray muons is applied using timing information from the 
muon system. 

Decays of r — > hadrons+^ T (called r^) are identified 
with a neural network [24 using as input variables (i) 
calorimeter clusters found with a cone algorithm of 1Z = 
0.3, (ii) energy in an annular cone 0.3 < 1Z < 0.5, 
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(Hi) electromagnetic (EM) calorimeter subclusters, (iv) 
the multiplicity of tracks with px > 1.5 GeV within 
1Z < 0.5 of the direction of the r lepton, and (v) con- 
sistency of the invariant mass of the hadron system with 
that of t decay. Three neural networks NN T are trained 
to identify tau decays corresponding to — > (j\), 
— > 7r^~7r®v (r 2 ), and -» 7r ± 7r ± 7r =F (7r°)i/ (r 3 ). In 
addition, a selection on the output of NN e i, a neural 
network trained to separate r 2 from electrons of similar 
signatures, is applied. The minimum transverse energy of 
the Th measured in the calorimeter, E T , is 12.5 GeV for 
T\ and t 2 and 15 GeV for r 3 . The sum of the transverse 
momenta of the r-associated tracks, p? k , is required to 
exceed 7, 5, and 10 GeV, for n, r 2 , and r 3 , respectively. 
In addition, at least one track with px > 7 GeV is re- 
quired for T3. Finally, only candidates with \rj Th \ < 1 
and Pt^ / E T > (0.65, 0.5, 0.5) for (n, r 2 , r 3 ) are selected. 

Jets are reconstructed from energies deposited in 
calorimeter towers using an iterative midpoint cone al- 
gorithm [35], with a cone radius 1Z = 0.5. Jet energies 
are calibrated to the particle-level jets using correction 
factors derived primarily from the transverse momentum 
balance in photon plus jets events Only jets with 

Py* > 15 GeV and \rj\ < 2.5 are considered in this anal- 
ysis and jets in the vicinity of a Th candidate (1Z < 0.5) 
are discarded. 

The Tf*T is calculated from the calorimeter energy, cor- 
rected for jet, EM, and r energy scales and for the trans- 
verse momentum of selected muons. 

Monte Carlo (MC) events for signal are simulated using 
MADGRAPH/madevent [27] and pythia [28] for parton- 
level generation and hadronization, respectively. We con- 
sider a range of scalar top quark mass values from 100 
to 200 GeV, generated in steps of 20 GeV. The range 
of probed sneutrino masses extends from 40 to 140 GeV 
in steps of 20 GeV. For each hypothesis, the MSSM pa- 
rameters are estimated from SUSPECT [29] and SDECAY 
[50] , The next-to-leading order (NLO) scalar top quark 
pair production cross section is calculated with prospino 
2.0 [3T], using CTEQ6.1M parton distribution functions 
(PDF) [32j [33] . The calculations are performed with the 
renormalization and factorization scales ^i r j equal to the 
scalar top quark mass roj , \ mj , and 2 rrij to estimate 
uncertainty on the nominal value through the impact of 
the two excursions. These uncertainties are combined 
quadratically with uncertainties on the PDF [33] [33] to 
provide a total theoretical uncertainty of 18% to 20% on 
the scalar top quark cross section. 

The kinematics for signal are determined both by 
njf and by the mass difference Am = m t ~ — mp. The 
Pt of the leptons and b quarks decrease on average for 
smaller values of Am, and Et is correlated with both 
mi and Am. We choose two signal points [rrif ,mp] = 
(180,60) GeV and (120,80) GeV, labeled "Signal A" and 
"Signal B" in the following, to illustrate the impact of 
the selection criteria for large mj and Am (Signal A) 



and for low mj and Am (Signal B). 

The dominant SM backgrounds to the pair produc- 
tion of scalar top quarks are from Zj^f*(-^ T + r~)+jets; 
z /l*(-^ ^ + M~)+j ct s; diboson production (WW, WZ, 
ZZ); it; VF+jets, and instrumental background from 
multijet (MJ) processes. All but the latter are estimated 
through MC simulations. Vector boson pair production 
is simulated with pythia, while the other backgrounds 
are simulated at the parton level using ALPGEN [33] and 
PYTHIA for hadronization and parton showering. 

Correction factors for MC estimated from data are ap- 
plied to lessen the impact of minor mismodeling of de- 
tector response. These corrections are related to the in- 
stantaneous luminosity, the position of the beam spot, 
identification efficiencies for \i and r, vector boson px, 
and jet, muon, and Th energy resolutions. 

The instrumental background originates either from 
incorrectly-identified, isolated muons (arising for exam- 
ple from semi-leptonic b decays) or from misidentificd 
Th (jets mimicking 77, signatures). This is estimated by 
changing the requirements on muon isolation and on 
the Th NN T outputs for each type of tau after the 
subtraction of the MC contributions corresponding to 
non-instrumental background. Normalization factors for 
these samples are estimated assuming that MJ processes 
have equal amounts of like-charge and opposite-charge 
[it events. 

The search of scalar top quark pairs proceeds in three 
steps: two event selections, labeled "Selection-1" and 
"Selection-2" below, and then a multivariate analysis. 

Selection-1 requires candidates to contain exactly one 
muon and 77, of opposite electric charge, and to have a 
minimum separation of 1Z([i,Th) > 0.5 between the two 
leptons. No specific requirement on jets is applied at this 
stage, but events having a jet in the vicinity of the muon 
(A7Z < 0.5) are rejected. Events with low values of 
the azimuthal angle difference between the leptons and 
Et, which are often due to issues with lepton reconstruc- 
tion, are removed by requiring A<p(fi, J$t) > 0.5 and 
Acj)(Th, St) > 0.5. At this stage of the analysis, 3387 
data events remain while 3453 ± 29 (stat) ± 440 (syst) 
events are expected from background. The main back- 
ground is from Z/j* (— > T + T~)+jets, M^+jets, and MJ 
events, as can be seen in Table [TJ Signal efficiencies do 
not exceed 4% for large Am, and are lower than 0.1% for 
Am < 20 GeV. 

Jets in scalar top quark pair events originate 
mainly from the hadronization of b quarks, whereas in 
Z/j*(— > T + T~)+jets and VF+jets backgrounds jets cor- 
respond predominantly to initial-state gluon radiation 
that provides a lower jet multiplicity. To maintain sen- 
sitivity to low Am signals, while rejecting a substantial 
part of the background, at least one jet is required in each 
event, which corresponds to Selection-2. Fig.[2]shows the 
jet multiplicity for the different Th decays. After this se- 
lection, 893 events remain, while a total background of 



() 



TABLE I: Numbers of events observed and expected from SM background processes and the two signal samples A and B at 
the Selection-1 stage for th (ti, T2, T3) and their sum. The uncertainties quoted in the Table are statistical. 



Process 


T\ 


T 2 


T3 


all r 


Z/7*(-> t+t~ )+jets 


162.6 


994.4 


352.3 


1509.4 


Z/7*(-> /i + M")+jcts 


38.7 


91.6 


48.3 


178.6 


diboson 


7.5 


40.5 


16.9 


65.0 


ti 


3.2 


27.3 


10.7 


41.3 


W+jets 


125.1 


631.1 


421.1 


1177.2 


Instrumental 


54.7 


233.1 


193.4 


481.3 


Background total 


392.0 ± 6. (J 


2018.1 ± 16.9 


1042.8 ± 12.6 


3452.8 ± 28.6 


Data 


388 


1937 


1062 


3387 


Wino scenario 










Signal A 


2.3 ±0.3 


17.4 ±0.8 


4.7±0.4 


24.4 ±1.0 


Signal B 


4.4 ±1.3 


21.4±3.0 


7.1 ±1.7 


32.9 ±3.7 


Higgsino scenario 










Signal A 


2.9 ±0.3 


20.5 ±0.9 


5.3±0.5 


28.6±1.1 


Signal B 


3.8±1.2 


20.5 ±2.9 


6.9±1.6 


31.1±3.6 



l in 3l (a) "tt 
2 10 !D0,L = 7.3fb- 1 □Z/ 7 * + jets 
w W+jets 
o 7^ ■ Inst. Bkg. 

□ Diboson 
• Data 

— Signal A 

— Signal B 




(b) 

D0, L = 7.3 fb 



23456789 10 
N(jets) 




□ Z/y*+jets 
HI W+jets 

Inst. Bkg. 

□ Diboson 
' Data 

— Signal A 

— Signal B 



1 23456789 10 
N(jets) 




FIG. 2: (Color online) Distributions of the number of jets after the preselection for (a) r type 1, (b) r type 2, (c) r type 3, 
and (d) their event sum. Signal A and Signal B correspond to the wino scenario. 



905 ± 10 (stat) ± 127 (syst) events is expected (see Table 




The separation of scalar top quark signal and back- 
ground is improved through an implementation of 
Boosted Decision Trees (BDT) [35] . A decision tree clas- 
sifies events on the basis of cumulative selection crite- 



ria that define disjoint subsets of events, each with a 
different signal purity. The decision tree is redefined 
iteratively by creating subsets of events called nodes. 
Each node is split into two subsets on the basis of 
the strongest discriminant for that sample. An im- 
purity measure i is estimated for each node from the 
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TABLE II: Numbers of events observed and expected from SM background processes and for the two signal samples A and B, 
after the final selection on JV(jets) > 0. The quoted uncertainties correspond to statistical sources. 



Process 


n 




T3 


all r 


JVQets) = 1 A(jets) = 2 AT (jets) > 2 


Z/7*(-> T + r~)+jets 


35.6 


226.2 


79.5 


341.3 


301.6 


34.6 


5.0 


Z/y*(-> M + /0+jets 


8.0 


20.7 


10.4 


39.1 


33.5 


4.7 


0.9 


diboson 


2.3 


12.5 


6.9 


21.7 


17.7 


3.3 


0.7 


ft 


3.0 


25.6 


10.1 


38.7 


8.1 


20.0 


10.6 


W+jets 


30.4 


166.5 


117.5 


314.4 


301.6 


36.7 


6.0 


Instrumental 


20.6 


55.0 


74.1 


149.8 


122.4 


20.1 


7.3 


Background total 


99.9 ± 2.3 


506.6 ± 5.4 


298.5 ± 4.9 


905.1 ± 9.6 


755.1 ±8.4 


119.4±2.0 


30.6 ±0.8 


Data 


90 


532 


271 


893 


738 


116 


39 


Wino scenario 
















Signal A 


1.9±0.3 


13.9±0.7 


3.7±0.4 


19.5 ±0.9 


10.4±0.6 


7.1 ±0.5 


2.0±0.3 


Signal B 


2.0±0.9 


10.6±2.1 


3.6 ±1.2 


16.3±2.6 


12.2 ±2.3 


3.1±1.1 


1.0±0.6 


Higgsino scenario 
















Signal A 


2.3±0.3 


16.2 ±0.8 


4.2±0.4 


22. 7± 1.0 


11.8±0.7 


8.6±0.6 


2.3±0.3 


Signal B 


1.8 ±0.8 


9.9±2.0 


3.1 ± 1.1 


14.8±2.4 


10.9 ±2.0 


2.9±1.0 


1.0±0.6 



weighted number of signal S and background B events 
in the node. For a given split, the decrease of impurity 
Ai = i(S, B) - i(S L , B L ) - i{S R , B R ), where L and R 
stand for left and right daughter nodes, is calculated. The 
best splitting gives the largest Ai. We use the Gini index 
[55] defined as SB/(S ± B) 2 as a measure of impurity. 
Terminal nodes are called leaves. Each leaf has a purity 
value defined by S/ (S + B). One of the main advantage 
of decision trees over analyses using simple requirements 
is that events that fail any individual selection criteria 
continue to be considered by the algorithm. 

The performance of the decision tree is improved by the 
boosting technique [36 . The basic principle is to create 
a tree, calculate an associated uncertainty, and create a 
new tree with a smaller uncertainty by re-weighting the 
misclassified events. We use adaptative boosting, known 
in the literature as AdaBoost [36] . The associated uncer- 
tainty e„ of a tree indexed by n is estimated as the frac- 
tion of misclassified events. The boosting weight of the 
n th tree is a n = /3 ln((l — e„)/e„) where f3 is an empir- 
ically determined parameter called boosting parameter. 
Misclassified events are given an additional multiplica- 
tive weight of e a " and the resulting new tree indexed by 
n + 1 is used to retrain the BDT and reduce the num- 
ber of misclassified events. This procedure is repeated N 
times, where N is the number of boosting cycles. For the 
scalar top quark search, the best signal to background 
separation occurs for ft = 0.5 and N = 40. 

To optimize the sensitivity of the analysis, three sub- 
samples are selected according to the jet multiplicity per 
event: iV(jets)= 1, AT(jets)= 2, N(jets)> 2. Since there 
are many SUSY mass points, and their characteristics 
differ significantly, the generated BDTs are trained and 
tested for each sample and each SUSY point using the 
implementation of the TMVA [37] library. The five most 
sensitive input variables for the BDTs trained with sam- 
ples of Signal A and Signal B are given in Table |III| To 
minimize bias, samples are split into three parts: 1/3 for 



training, 1/3 for testing and 1/3 for analysis. The distri- 
butions of the BDT outputs trained with Signals A and 
B are shown in Figs. [3] and [4] for the wino and higgsino 
scenarios, respectively. 

The predicted numbers of background and signal 
events depend on measurements and parametrizations 
that have non-negligible systematic uncertainties, which 
can either affect exclusively the normalizations of back- 
grounds or the signal efficiency, or modify also the differ- 
ential distribution of the BDT discriminant. The main 
sources involve muon identification and reconstruction ef- 
ficiencies (2%), r h identification and reconstruction (10%, 
4%, and 5% for n, T2 and T3, respectively), trigger (5%), 
luminosity (6.1%) [38], jet energy calibration (3.2% for 
the background, 1.5% to 2.6% for signal), jet identifica- 
tion efficiency and energy resolution (5% for the back- 
ground, 1% to 6% for signal). Systematic uncertainties 
related to reconstructed objects are estimated by chang- 
ing each quantity by one standard deviation (s.d), and 
gauging the impact on the final measurement. Additional 
uncertainties arise from the choice of PDF, which affects 
the cross sections for the background components (6.3% 
for Z/~/*{-^ r+T-)±jcts and Z/ 7 *(^ /i+>-)+jets, 15% 
for W+]ets and 10% for ft; 5.6%, 8.1%, and 5.5% for 
WW, WZ, and ZZ, respectively) and for signal (18% 
to 20%). To estimate the systematic uncertainties re- 
lated to instrumental background, a MC scalar top quark 
signal is added as a background contribution during the 
stages that consider distribution and normalization of the 
MJ background. Each scalar top quark signal is consid- 
ered, and the largest relative changes in the distribution 
and the normalization of the instrumental background 
are thereby estimated. The resulting uncertainties cor- 
respond to 5% and 10% for the normalization and the 
functional dependence, respectively. 

There is no significant excess of events observed above 
the predicted background, and we therefore combine the 
numbers of expected signal and background events, with 



TABLE III: Listing of the five most sensitive input variables used for each BDT training and testing. The significance of 
Bt, Sig(^r), is defined as the likelihood that the Bt in an event is consistent with a fluctuation of the resolution on the 
Pt measurements on the selected leptons and jets. St is the sum of the lepton pr and of the Bt- The transverse mass Mt is 
defined as Mt(A, B) = y/2p^p^, (1 — cos A(j>(A, B)). Ht is equal to the scalar sum of the Et of the jets. 





Wino scenario 






Higgsino scenario 


JV(jete)= 1 


N (jets) = 2 


JV(jets)> 3 


A(jets) = 1 


iV(jets)= 2 


JV(jets)> 3 


S T 




A0 mln (jet,^ T ) 


A<T ln (jet,9 T ) 


Mass(^i, r) 


A0 m ' n (jet,^ T ) 


^(leading jet) 


^(leading jet) 


Sig(B T ) 


^(leading jet) 


^(leading jet) 


Ht 


Signal A Aft max O,jet) 


Ar ax (rjet) 


Aft max (Tjet) 


Aft max (r,jet) 


Aft max (A*Jet) 


Aft max (r,jet) 


M T (fJ,,B T ) 


A<j>(nr,B T ) 




Sig(B T ) 


A<j>(n + r, B T ) 


Mass(^i, r) 


Mass(lept Jet) 


A7l max ( M Jet) 


Pt (leading jet) 


Mass (lept Jet) 


M t (t, B t ) 


Mass(/ijet) 


Sig(9 T ) 


A0 mln (jet,$ T ) 


A^Qet,B T ) 




A<r in (jet, B T ) 


A0(/i,r) 


T)(t) 


^(leading jet) 


Sig(B T ) 


St 


^(leading jet) 


A<^(leading jet,^T) 


Signal B S T 


Aft mi >,jet) 


A<t>{pL + T, B T ) 


n(r) 


Aft min (r,jet) 


Aft max (rjet) 




A0(next-to-lcading ]et,BT) 


A (^(leading ]et,Br) 


A<£0,t) 


A^(lcading jot, Bt) 


A4>{n,B T ) 


Mass(lept Jet) 


A7l min ( M Jet) 


Mass(/ijet) 


Mass (lept Jet) 


A0( M + r, B T ) 


A(t>(r, B T ) 



(a) 

D0, L = 7.3 fb" 1 
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□ z/y*+jets 
W-fjets 
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■ Data 

— Signal A 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
BDT output 



:7.3 fb" 1 
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□ z/y*+jets 
Hw+jets 
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BDT output 
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LU 



(e) 

D0, L = 7.3 fb" 1 



■ tt 

□ z/y*+jets 
W+jets 
Inst. Bkg. 

□ D Iboson 
• Data 

Signal A 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
BDT output 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 lu 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

BDT output BDT output BDT output 



FIG. 3: (Color online) Distributions of the BDT output discriminants, in the wino scenario, for the sample with A r (jets)= 1, 
(a) for Signal A, (b) for Signal B; iV(jets)= 2, (c) for Signal A, (d) for Signal B; A r (jets)> 2, (e) for Signal A, (f) for Signal B. 



their corresponding uncertainties, and the numbers of 
events observed in data obtained from the BDT outputs 
for each SUSY point, to calculate upper limits on the 
cross sections for signal at the 95% CL using the modi- 
fied frequentist approach [3S] . The bins of the BDT out- 
puts corresponding to iV(jets)= 1, 2, and > 2, are treated 
as separate channels, and their likelihoods are combined 



taking into account the correlations of both systematic 
uncertainties affecting exclusively the normalization of 
backgrounds and signal efficiencies and also of those that 
change the distribution of the BDT discriminant. The 
limits are calculated using the confidence level CLg — 
CLs+b/CLb where CLs+b and CLb are the confidence 
levels for the signal+background and background-only 
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W+jets 
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BDT output 



10 



(e) 

aLD0, L= 7.3 fb" 



10 2 : 



■ ft 

HI Z/y*+jets 
H W+jets 

Inst. Bkg. 

Dlboson 
• Data 
— Signal A 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
BDT output 
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BDT output BDT output BDT output 



FIG. 4: (Color online) Distributions of the BDT output discriminants, in the higgsino scenario, for the sample with 7V(jets) = 
1, (a) for Signal A, (b) for Signal B; A r (jets)= 2, (c) for Signal A, (d) for Signal B; 2V(jets)> 2, (e) for Signal A, (f) for Signal 
B. 



hypotheses, respectively [35]. Exclusion regions are given 
on Figs. [5] and § as a function of the scalar top quark and 
sneutrino masses. These results are obtained under the 
assumption B(t\ -> b\w) = B(ti -> brv) = 1/3 (Fig. [H] 
wino scenario) and B(ti — > bfiu) = 0.1, B{t\ ~ > brv) 
— 0.8 (Fig. [6j higgsino scenario). For larger mass differ- 
ences between the scalar top quark and the sneutrino, a 
scalar top quark mass lower than 200 GeV is excluded. 
The search is sensitive to a possible signal in the mass 
region up to Am = 60 GeV for = 140 GeV, with 
the observed limit being within one standard deviation 
of the expected limit. 

In summary, a search for scalar top quark pair pro- 
duction in pp collisions at \fs = 1.96 TeV has been per- 
formed in a dataset corresponding to an integrated lu- 
minosity of 7.3 fb -1 . Events containing one muon, one 
r decaying hadronically, at least one jet, and missing 
transverse energy have been considered in this analysis. 
No evidence is found for the production of the lightest 
scalar top quark, and 95% CL exclusion limits are set in 
the plane [mj ,mp]. The largest scalar top quark mass 
excluded is 200 GeV for a sneutrino mass of 45 GeV, 
and the largest sneutrino mass excluded is 85 GeV for a 



scalar top quark mass of 160 GeV. This is the first Teva- 
tron limit obtained from a study of final states containing 
r leptons from t\t\— > bbfirB^T- 
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